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Abstract 

We study theoretically the spatiotemporal response of a lipid membrane submitted to a local chemical change of its environment, 
taking into account the time-dependent profile of the reagent concentration due to diffusion in the solution above the membrane. 
We show that the effect of the evolution of the reagent concentration profile becomes negligible after some time. It then becomes 
possible to extract interesting properties of the membrane response to the chemical modification. We find that a local density 
asymmetry between the two monolayers relaxes by spreading diffusively in the whole membrane. This behavior is driven by 
intermonolayer friction. Moreover, we show how the ratio of the spontaneous curvature change to the equilibrium density change 
induced by the chemical modification can be extracted from the dynamics of the local membrane deformation. Such information 
cannot be obtained by analyzing the equilibrium vesicle shapes that exist in different membrane environments in light of the area- 
difference elasticity model. 
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1. Introduction 

During cell life, membranes are submitted to an inhomo- 
geneous and variable environment. Local inhomogeneities can 
be strongly related to biological processes, which has led to 
experiments investigating the effect of local modifications on 
biomimetic membranes jlj]. For instance, in the inner mem- 
brane of mitochondria, the enzymes that use the local pH dif- 
ference across the membrane to synthesize adenosine triphos- 
phate, the cell's fuel, are located in membrane invaginations 
called cristae [21]. Experiments on model lipid membranes have 
shown that a local pH change can induce a local dynamical 
membrane deformation 0, El Q 01, and in particular the 
formation of cristae-like invaginations [3]: membrane shape is 
tightly coupled to local pH inhomogeneities. Other concentra- 
tion inhomogeneities in the environment of a cell have a crucial 
biological role, for instance in chemotaxis or in paracrine sig- 
naling. It is therefore of great interest to study the response of a 
biological membrane to a local modification of its environment. 

Motivated by experiments conducted on biomimetic mem- 
branes, we have developed a theoretical description of the dy- 
namics of a lipid bilayer membrane submitted to a local con- 
centration increase of a substance that reacts reversibly and in- 
stantaneously with the membrane lipid headgroups. We focus 
on the regime of small deformations, and we treat linear mem- 
brane dynamics in the spirit of Ref. |8|]. While our first works 
focused on the simple case of a constant modification of the 
membrane involving only one wavelength JiHH, we have re- 
cently extended our theoretical description in order to take into 
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account the spatiotemporal profile of the fraction of chemically 
modified lipids resulting from the local reagent concentration 
increase |2[. This profile is determined by the diffusion of 
the reagent in the solution that surrounds the membrane. In 
Ref. j7|], we compared the predictions of this theoretical de- 
scription to experimental measurements of the deformation of 
the membrane of giant unilamellar vesicles caused by local mi- 
croinjection of a basic solution, and we obtained good agree- 
ment between theory and experiment. 

In the present article, we pursue the theoretical investigation 
of the effect of a local chemical modification on a lipid mem- 
brane. In general, the dynamics that results from a local reagent 
concentration increase is quite complex, as it involves the evo- 
lution of the reagent concentration profile simultaneously as the 
response of the membrane. This is the case in the experimen- 
tal data analyzed in Ref. [7], which corresponds to microin- 
jection steps lasting a few seconds. Here, we show that the 
effect of the evolution of the reagent concentration profile on 
the membrane dynamics becomes negligible some time after 
the beginning of the reagent concentration increase, after what 
the dynamics corresponds to the response of the membrane to 
a chemical modification imposed instantaneously. We find that 
studying this regime enables to extract interesting properties of 
the membrane response. 

The article is organized as follows. First, in Sec. |2l we re- 
view the linear dynamics of a membrane submitted to a local 
chemical modification. Then, in Sec. [3] we study separately 
the dynamics associated with each of the two effects that can 
arise from a chemical modification, namely a spontaneous cur- 
vature change and an equilibrium density change of the external 
monolayer. We find that a local asymmetric density perturba- 
tion between the two monolayers of the membrane relaxes by 
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spreading diffusively in the whole membrane. Intermonolayer 
friction plays a crucial part in this behavior. Subsequently, in 
Sec. |U we treat the general case where both effects are present, 
and we show how the ratio of the spontaneous curvature change 
to the equilibrium density change induced by the local chemi- 
cal modification can be extracted from the dynamics. This ratio 
cannot be deduced from the study of global modifications of 
vesicle equilibrium shapes in light of the area-difference elas- 
ticity model [9]. Finally, Sec.[5]is a conclusion. 

2. Dynamics of a chemically modified membrane 

For the article to be self-contained, the present section re- 
views the linear dynamics of a membrane submitted to a local 
chemical modification, starting from first principles. The main 
points of this description were presented in Ref. [7]. In that ar- 
ticle, we compared theoretical predictions to experimental mea- 
surements of the deformation of a membrane submitted to a lo- 
cal and brief pH increase. Here, our aim will be to go further in 
the analysis of our theoretical description in order to understand 
the fundamental properties of the response of a membrane to a 
local chemical modification. 

2.1. Monolayer free energy 

Our description of the bilayer membrane is based on a lo- 
cal version of the area-difference elasticity membrane model |Hl 

Ml 

1 ill . We focus on small deformations of an infinite flat mem- 
brane, and we denote the upper monolayer by + and the lower 
one by -. 

In the absence of a chemical modification, the local state 
of each monolayer is described by two variables: the local to- 
tal curvature c defined on the membrane midlayer, which is 
common to both monolayers, and the local scaled density r ± = 
D± - Po)/po, defined on the midlayer of the membrane, po be- 
ing a reference density. The sign convention for the curvature 
is chosen in such a way that a spherical vesicle has c < 0. The 
free energy /* per unit area in monolayer + reads 111 ill : 



cr K 2 KC k 2 
f — — — I- —c + —c + -(r ± ± ec) 
2 4 2 2 



(1) 



where cr represents the tension of the bilayer and k its bend- 
ing modulus, while k is the stretching modulus of a monolayer, 
and e denotes the distance between the neutral surface 1 1 211 of 
a monolayer and the midsurface of the bilayer. As we assume 
that the two monolayers of the membrane are identical before 
the chemical modification, these constants are the same for both 
monolayers. The spontaneous curvatures of the two monolayers 
have the same absolute value cq and opposite signs, since their 
lipids are oriented in opposite directions. The expression for f* 
in Eq. ([TJ corresponds to a general second-order expansion in 
the small dimensionless local variables r ± and ec, around the 
reference state which corresponds to a flat membrane with uni- 
form density p ± - po. It is valid for small deformations around 
this reference state: = 0(e) and ec = 0(e), where e is a 
small dimensionless parameter used for bookkeeping purposes, 
which characterizes the amplitude of the small deformations of 



the membrane around the reference state. Mathematically, e is 
considered infinitesimal. Note that in general, both c and r*, 
which describe local small deformations around the reference 
state, are functions of time and of position on the membrane. 

Let us now focus on the way the membrane free energy is 
affected by the local chemical modification. We consider that 
the reagent source, which corresponds to the micropipette tip 
in an experiment, is localized in the water above the mem- 
brane. Besides, membrane permeation and flip-flop are ne- 
glected given their long timescales. Hence, the chemical modi- 
fication only affects the upper monolayer, i.e., monolayer +, and 
not the lower one. Let us denote by <p the local mass fraction of 
the lipids of the upper monolayer that are chemically modified: 
<p depends on time and position since it arises from the local 
chemical modification. We assume that the reagent concentra- 
tion is small enough for <p to remain small at every time and 
position on the membrane, and we characterize this smallness 
through <p = 0(e). In order to describe the chemically modified 
membrane, we have to include the third small variable <p in our 
second-order expansion of / + . We obtain 111 111 : 



k k , 2 

+ 2 ( c o + co4>) c + - (r + ec) , 



(2) 



where the constants a~\, o~2, and cq describe the response of 
the membrane to the chemical modification. These constants 
depend on the reagent that is injected. Their physical mean- 
ing will be explained in the next paragraph. Besides, the non- 
analytical mixing entropy term cr(l + r + ) <pln<p has been added 
to our second-order expansion ll III . Note that we assume that 



the three small dimensionless local variables <p, r* and ec are 
of the same order. In fact, in the present work, the deformation 
of the membrane and the density variation are caused by the lo- 
cal chemical modification, i.e., they are a response to <p, which 
justifies that ec and r ± are of the same order as (p. We refer the 
reader to Ref. fllll for more details on the derivation of Eqs. ([TJ 
and ©. 

The effect of the chemical modification (i.e., of (p) on the 
upper monolayer is twofold. First, the scaled equilibrium den- 
sity on the neutral surface of the upper monolayer is changed by 
the amount cr\<pjk to first order. Second, the spontaneous curva- 
ture of the upper monolayer is changed by the amount -co4> to 
first order, with cq = Co + 2cr 1 e//e. These results are obtained by 
minimization of the free energy of a homogeneous monolayer 



with constant mass (see Appendix A i. Hence, the constants cr\ 
and co describe the linear response of the monolayer equilib- 
rium density and of its spontaneous curvature, respectively, to 
the chemical modification. This explains the physical meaning 
of the constants cr\ and cq in Eq. (0. Note that o~2 corresponds 
to the quadratic response of the membrane to the chemical mod- 
ification, but it will not have any relevant effect in the following. 

2.2. Dynamical equations 

The elastic force densities in a monolayer described by the 



free-energy densities in Eqs. (QJ[2| have been derived in Ref. Ill 111 
to first order in e, using the principle of virtual work. As we 
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focus on small deformations of an infinite fiat membrane, it 
is convenient to describe it in the Monge gauge by the height 
z = h{r), r e R 2 , of its midlayer with respect to the reference 
plane z — 0. Then, ec = eV 2 h to second order. Such a descrip- 
tion is adapted to practical cases where the distance between 
the reagent source and the membrane is much smaller than the 
vesicle radius. The force per unit area of the reference plane, 
which we call "force density", then reads to first order in e 

pf = -k V (r + + eV 2 h -j^J, (3) 
p- = -kVir -el 2 h) , (4) 
p z = cr V 2 /z - kl A h - ke l 2 r a - - cr iej V 2 , (5) 

where pf is the tangential component of the force density in 
monolayer "±", while p z = pi + p~ is the total normal force 
density in the membrane. In these formulas, we have introduced 
the antisymmetric scaled density r a — r + - r~ , and the constant 
k — k + 2ke 2 . These expressions show that both the equilibrium 
density change and the spontaneous curvature change (i.e., both 
0"i and Co) can yield a normal force density, and thus a defor- 
mation of the membrane, while only the equilibrium density 
change can yield a tangential force density and induce tangen- 
tial lipid flow. An illustration of the role of these force densities 
in the membrane response to a local chemical modification is 
provided in Fig. Q] 
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Figure 1 : Qualitative representation of the early stages of the membrane re- 
sponse to an instantaneous chemical modification affecting the upper mono- 
layer (monolayer +). (a): Membrane at equilibrium in the reference state (flat 
shape, uniform lipid densities in both monolayers), (b): Membrane just after 
the instantaneous local chemical modification: the membrane has not deformed 
yet at this stage, and the densities are still uniform. Some lipids of monolayer 
+ are modified: they are represented in dark blue (<j> represents their local mass 
fraction; for clarity <p is locally large in the figure, but in our work, cj> = 0(e)). 
Here we assume that the result of the chemical modification is to increase the 
effective size of their headgroups. Hence, the equilibrium density in monolayer 
+ is locally decreased (o"i < 0), while the spontaneous curvature is increased 
in absolute value (cq > 0). This results into a normal force density p z > at the 
center of the modified zone of the membrane (see Eq. (5))> and into a tangential 
force density in monolayer +, oriented toward the exterior of the modified zone 
(see Eq. (5J). Both are indicated on the figure by arrows, (c) Due to the normal 
force density, the membrane starts deforming upwards. The maximum height 
H and the full width at half maximum W of the deformation are indicated. Note 
that tangential flows also appear due to tangential force densities. 



to describe the dynamics of the membrane to first order in the 
spirit of Ref . {§] . Let us outline the derivation of these equa- 
tions JHQl, which is described in more detail in Appendix B 



We use a normal force balance for the bilayer (see Appendix B.2 1 
which involves p z and the normal viscous stress exerted by the 
surrounding fluid, which is derived in Appendix B.l Besides, 
as each monolayer is a two-dimensional fluid, we write down 
generalized Stokes equations (see | Appendix B.2[ ), which in- 
volve pf, the two-dimensional viscous stress associated with 
the lipid flow, the tangential viscous stress exerted by the sur- 



rounding fluid (see Appendix B.l i, and the intermonolayer fric- 
tion +b (v + - v~), where v* is the velocity in monolayer + and 
b is the intermonolayer friction coefficient H, 13]. This inter- 
monolayer friction is a tangential force that comes into play 
when one monolayer slides relatively to the other on^. Finally, 
we use mass conservation in each monolayer. These dynami- 
cal equations are best expressed using two-dimensional Fourier 
transforms of the various fields involved, denoted with hats: for 
any field / which depends on r and on time t, f is such that 



f(q,f)= drf(r,f)e- i " r . 



(6) 



Combining all the above-mentioned equations yields a system 
of first-order linear differential equations on the two-dimensional 
variable X = (qh, r a ): 

^(q,t) + M(q)X{q,t) = Y(q,t), (7) 
where we have introduced the matrix which describes the dy- 



namical response of the membrane [8l |5| 



M(q) = 



cr q + kq 3 


keq 


4i] 


4i] 


keq 3 


kq 2 


b 


~2b 



(8) 



Here, we have assumed that r]2q 2 <Z b and rjq <K b. This is true 
for all the wavevectors with significant weight in 0, if the modi- 
fied lipid mass fraction (f> has a smooth profile with a character- 
istic width larger than 1 /im. Indeed, 77 = 10~ 3 J s/m 3 for water, 
and typically 772 = lO^Js/m 2 and b = 10 9 Js/m 4 0(3. 
Besides, the forcing term in Eq. (0 reads J^: 

-,2 



Y(q,t) = 



KCoq 



o~\q 

2b 



<P(q,t) 



(9) 



Eqs. (|7][8]) show that the membrane deformation is coupled 
to the antisymmetric density: the symmetry breaking between 
the monolayers causes the deformation of the membrane. Here, 
the symmetry breaking is caused by the chemical modification 
of certain membrane lipids in the external monolayer, i.e., to 
the presence of (p. And indeed, Eq. (0 shows that the forcing 
term in Eq. (0 is proportional to (j>(q, f). 



Using the elastic force densities in Eqs. OHSJ), it is possible tangential 



2 Note that the two monolayers are supposed to remain in contact with one 
another at every time and point, due to hydrophobic forces: the sliding is only 



3 



2.3. Profile of the fraction of chemically modified lipids 

Before solving Eq. (0, we need to determine cf>(q, t), which 
is involved in the forcing term Y(q, t). The profile <p{r, t) of 
the mass fraction of chemically modified lipids in the exter- 
nal monolayer arises from the local reagent concentration in- 
crease. We focus on reagents that react reversibly with the 
membrane lipid headgroups. Besides, we assume that the re- 
action between the lipids and the reagent is diffusion-controlled 
(see, e.g., Ref. flaO . In other words, the molecular reaction 
timescales are very small compared to the diffusion timescales. 
For such a reversible diffusion-controlled chemical reaction, 
<f>(r, t) is instantaneously determined by the local reagent con- 
centration on the membrane, which results from the diffusion 
of the reagent in the fluid above the membrane. Note that these 
hypotheses are verified in the experiments analyzed in Ref. |l7J], 
where the reagent is sodium hydroxide. 

We consider that the reagent source is localized in (r, z) = 
(0, zo > 0), which would represent the position of the micropipette 
tip in an experiment (see Ref. |01). The cylindrical symme- 
try of the problem then implies that the fields involved in our 
description only depend on r = \r\. We focus on the regime 
of small deformations h(r, t) <K zo, and we work at first order 
in h(r, t)/zo- Besides, we study the linear regime where <p{r, t) 
is proportional to the reagent concentration on the membrane: 
denoting by C(r,z,t) the reagent concentration field, we have 
<p(r, t) oc C(r, h(r, /), f). To first order in the membrane deforma- 
tion, this can be simplified into <p(r, t) oc C(r, 0, t). 

The field C is determined by the diffusion of the reagent 
from the local source in the fluid above the membrane. In mi- 
croinjection experiments, there is also a convective transport of 
the reagent due to the injection, but the Peclet number remains 
so small that diffusion dominates. Besides, since the membrane 
is a surface and <p <k 1, we can neglect the number of reagent 
molecules that react with the membrane when calculating C. 
Hence, C can be obtained by solving the diffusion equation 



d,C 



D1 2 C 



(10) 



where D is the diffusion coefficient of the reagent in the fluid 
above the membrane, and the source term reads 



S(r,z,i) = S S(r)6(z-ZoMi), 



(11) 



where is Heaviside's function. This corresponds to a con- 
stant injection flow from the source starting at t = 0. In ad- 
dition, the membrane imposes a Neumann boundary condition 
d z C (r, h(r, t), t) = 0. This relation, which corresponds to a van- 
ishing flux across the membrane, can be simplified to first order 
into 

d z C(r,0,t) = 0. (12) 
The solution to this diffusion problem reads 

pt p p+oo 

C (r, z,t)= I df I dr' I dz' S (r',z, ?') x 
Jo Js? Jo 

G(\r-r'\,z,z',t-t') 

= S f dt'G(r,z,z Q ,t-t') , (13) 
Jo 



where the causal Green function G of our diffusion problem can 
be expressed using the method of images ifnll : 

G (r, z, z', t) = G ra (r, z - z' , t) + G m {r, z + z' , t) , (14) 

where we have introduced the infinite-volume causal Green func- 
tion of the diffusion equation 



Geo (r, z, t) = 



6(t) 



{AirDtf 12 



exp - 



r 2 + z 2 
\Dt 



(15) 



Combining Eqs. ( TT3l fT4l [TBT l provides an analytical expression 
for C(r, z, t), and for C(q, z, t). Since (p(q, t) oc C (q, 0, t), we thus 
obtain an analytical expression for 0, which reads 

^,0ac«f(,VS--5 = )25^-^^, (16) 

\ i^Dtl qzo qzo 

where erf denotes the error function. 
As t — > oo, (j>(q, t) converges towards 



<f>s(q) ^ 



-qz u 



qzo 



(17) 



2.4. Resolution of the dynamical equations 

The membrane deformation h is given by the solution to 
Eq. (|7]i with expressed in Eq. ( flGb and with the initial con- 
dition X(q, t = 0) = (0, 0), corresponding to a non-perturbed 
membrane. Hence, we have 

qh(q,t) = ds [vie- ri 'A(s) + v 2 e~ r2 'fi(s)] . (18) 



Here, we have introduced the eigenvalues j\ and 72 of the ma- 
trix M(q) (see Eq. (O), the associated eigenvectors V\ = (vi, w{) 
and V2 = (V2, W2), and the solutions A and B of the linear system 



ViA(t)e- y " + V 2 B(t)e 



-72' 



Y(q,t). 



(19) 



The integral in Eq. ( fT8l ) and the inverse Fourier transform of h 
can be calculated numerically. This gives the spatiotemporal 
evolution of the membrane deformation h. 

2.5. Dimensionless form of the description 

Both for physical understanding and for computational effi- 
ciency, it is interesting to put our description in dimensionless 
form. Using the Buckingham Pi theorem and choosing zo as the 
distance unit, bz\jk as the time unit and a-o^z^/lc 2 as the mass 
unit, the ten parameters in our equations yield seven dimension- 
less numbers: 



zo V °"o ' zo ' ° k 



Db o- bzo 

A = — , p = 

k r\k 



£1 = -7- , G = cozo . (20) 
k 



The parameters L\ and L 2 compare the characteristic length- 
scales ^k/ctq and e of the membrane to zo, while So is a di- 
mensionless version of the membrane tension cr . Besides, the 
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parameter A quantifies the importance of the reagent diffusion 
on the membrane dynamics. We will discuss briefly the effect 
of varying A in the following. 

In the case of a lipid membrane in water, the only dimen- 
sional parameters that can span various orders of magnitudes 
are cr and zp- For cr > 10~ 8 N/m, i.e., for realistic membrane 
tensions 11181 Il9ll . and for zo > 5 /mi, we will see in the fol- 
lowing that the only relevant dimensionless parameters in the 
dynamics of the membrane are /! and a = -Y.qL 2 G /(L2E1). The 
first one, //, then corresponds to the ratio of the two eigenvalues 
of M(q) for q = l/(2zo) (see Sec. 13.1.21 ): hence, it is a crucial el- 
ement of the membrane dynamical response. The second one, 
a, quantifies the relative weight of the spontaneous curvature 
change and of the equilibrium density change of the external 
monolayer due to the chemical modification (see SecHJi. Note 
indeed that a oc G/Ei, and that G and Si describe the two ef- 
fects of the reagent on the membrane, through co and o~\ . The 
effect of varying p. and a will thus be discussed in the following. 

For our numerical calculations, we initially take 



10" 1 , L 2 



1(T 4 , E 



10 



-6 



A = 21.25,//= 10. 



(21) 



These values correspond to the injection of NaOH in water (D = 
2125 /mi 2 /s lEoll ) from a source at zo = 10/mi above a floppy 
membrane with typical constants <x = 10~ 7 N/m (see Ref. 11911 
and references therein), k = 10~ I9 J, k — 0.1 N/m and e — 
1 nm 1 12]. The values of the intermonolayer friction coefficient 
b and of the water viscosity r] are those given above. Note that 
the time unit bz^/k is then equal to one second. 

3. Extreme cases 

In this Section, we will study the two extreme cases G — 
and Si = 0. In other words, we will study separately the dy- 
namics associated with an equilibrium density change and with 
a spontaneous curvature change of the external monolayer. We 
will see that these two manifestations of the chemical modifica- 
tion of the membrane due to the reagent concentration increase 
yield different spatiotemporal evolutions of the membrane de- 
formation. The general case where both effects are present will 
then be discussed in Sec.|4] 

A priori, the dynamics of the membrane deformation is quite 
complex, as it involves the evolution of the reagent concentra- 
tion profile, due to diffusion, simultaneously as the response of 
the membrane. In order to investigate the effect of the reagent 
diffusion on the dynamics of the deformation, and to see when 
the effect of the evolution of the reagent concentration profile 
becomes negligible, we will compare the two following cases: 

(i) the realistic case where <p(q, t) is given by Eq. ([Tol l. 

(ii) the theoretical case where the stationary modification <p s in 
Eq. (fT7T > is imposed instantaneously: cf>{q, t) = (f> s (q)6(f). Note 
that this case corresponds to the limit of very large A. 

3.1. Equilibrium density change only 

Let us first focus on the case where G — 0, in which only the 
equilibrium density of the upper monolayer is affected by the 



chemical modification. Fig.|2]shows the evolution with dimen- 
sionless time r = kt / (bzfy of the membrane deformation height 
H(j) = h(0, r) in front of the source, and of the full width at 
half-maximum W(t) of the deformation (see Fig.QJc) for an il- 
lustration of the definitions of H and W). The cases (i) and (ii) 
introduced above are presented. On these graphs, we also show 
<E>(t) = 0(0, t)/</>. s (0) and the full width at half-maximum of 
4> in case (i). 

Fig.|2]shows that the membrane undergoes a transient defor- 
mation that relaxes to zero while getting broader and broader. 
The local asymmetric density perturbation relaxes by spreading 
in the whole membrane. However, this process is slowed down 
by intermonolayer friction, so the density asymmetry is tran- 
siently solved by a deformation of the membrane j2ll l5ll. This 
deformation is downwards if Ei > 0, and upwards if E] < 0. 

3.1.1. Effect of reagent diffusion 

In the realistic case (i), at short times, the membrane dy- 
namics is governed by the evolution of <p due to the reagent 
diffusion. Conversely, at long times, once </> is close enough to 
its steady-state profile <f> s , the dynamics of the membrane defor- 
mation is similar in case (i) and in case (ii). This can be seen 
in Fig. [2] first, H{\) and W(i) follow <t> and W$, and then, they 
have an evolution very similar to those of //(ii) and W(ii). Thus, 
after some time, the effect of reagent diffusion on the dynam- 
ics of the membrane deformation becomes negligible, and the 
dynamics of the realistic case (i) can be well approximated by 
that of the simpler case (ii), which corresponds to the response 
of the membrane to an instantaneously imposed modification. 

The transition time t c between the diffusion-dominated regime 
and the membrane-response-dominated regime is determined 
by the convergence of <p to <p s . As the reagent takes a dimen- 
sionless time 1/A to diffuse from the source to the membrane, 
we expect t c oc 1/A. We studied the dynamics of the defor- 
mation height for A € [10~ 3 , 10 3 ] by integrating Eq. © nu- 
merically, and we found that this law is very well verified. For 
our standard value of A (see Eq. l2TT i. used in Fig. [2] we have 
t c « 0.5 ~ 10/A. 

3.1.2. Analytical insight 

Let us discuss analytically the simple case (ii). The long- 
time behaviors obtained in this case are especially interesting, 
since they also apply to the realistic case (i). Let us focus on 
cr > 10~ 8 N/m, i.e., on realistic membrane tensions, and let 
us keep the standard values of the other parameters involved 
in M(q) (see below Eq. (OTl i). as these parameters cannot vary 
significantly for a membrane in water. The eigenvalues of M(q), 
which correspond to the two independent relaxation rates of a 
perturbation of the membrane with wavelength 2n/q, can then 
be approximated by 



yi 



kq 2 

' ~2b 

crpq 

Ar) 



(22) 
(23) 



for all wavevectors with significant weight in (f> s if zq > 5 /urn. 
Within this approximation, for q = l/(2zo), we have fi = 72/71- 
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Figure 2: Dynamics of the membrane deformation in the extreme case G = 0, 
where only the equilibrium density is changed. The values taken for the other 
dimensionless numbers are those in Eq. 12H . Both the realistic case (i) where 
reagent diffusion is accounted for, and the simpler case (ii) where the chemical 
modification is imposed instantaneously, are considered, (a) Logarithmic plot 
of the height of the membrane deformation H and of the fraction <I> of modified 
lipids in front of the reagent source versus dimensionless time r. Both in case 

(i) and in case (ii), H is plotted in units of the extremal value it attains in case 

(ii) . (b) Logarithmic plot of the width W of the membrane deformation and of 
the width of the fraction of modified lipids versus r. Both W and W<j are 
plotted in units of zo- It can be seen on graphs (a) and (b) that cases (i) and (ii) 
yield similar dynamics for r 2> r c x 0.5. The thin red (gray) lines correspond 
to the analytical laws mentioned in the text. 



This leads to a simple expression of the solution of the dynam- 
ical equation Eq. (0: 



e -p/Ji"/4 _ e -p 2 r/2 
Hp, T) oc $ s (p) , 

P ~2p 



(24) 



with p = qzo- This expression shows that the only parameter 
that is relevant in the dynamics is p. Note that, in the realistic 
case (i), the value of A is relevant too. 

In the long-time limit, as the deformation spreads, we have 
2p <sc p for all the wavevectors with significant weight in h(p, t). 
In this case, calculating the inverse Fourier transform of Eq. ( 1241 
for r — yields 

H-l/V?, (25) 
and the numerical results for W are in excellent agreement with 



W ~ V7, 



(26) 



as can be seen in Fig.[2|b). This law can also be obtained ana- 
lytically if <f> s is replaced by a Gaussian. 

In the short-time limit, Eq. (f24b yields H ~ r and W — > 



time limit and for the much more interesting long-time limit, 
are plotted in red (gray) lines on Fig. [2] 

The transition between the two asymptotic regimes is de- 
termined by the value of p, which is the only parameter that 
controls the dynamics in case (ii). In particular, it is possi- 
ble to show that the long-term scaling laws are valid for r » 
max (1, l/p 2 ). We studied the dynamics of the deformation 
height for fi e [1CT 3 , 10 3 ] by integrating Eq. (0 numerically, 
and we found that this law is very well verified. 

3.1.3. Diffusive spreading of an antisymmetric density pertur- 
bation 

The long-term scaling W ~ -\Jt shows that the local anti- 
symmetric density perturbation spreads diffusively. Let us em- 
phasize that this diffusive behavior is not related to the diffusion 
of the reagent in the solution above the membrane, as in the 
long-term limit, <p has reached its steady-state profile <p s . The 
fact that the long-term scaling W ~ yfr holds in case (ii), where 
the profile <p s is established instantaneously, as well as in case 
(i) (see Fig. HJb)), illustrates that this scaling law is not related 
to the reagent diffusion. 

The long-term diffusive spreading of the antisymmetric den- 
sity perturbation can be understood as follows. At long times, 
the difference between the Stokes equations for each monolayer 
Eq. JB. 12b and Eq. ( IB. 1 3b can be approximated by 

kv(r a -^-c^ + 2b{v + -v-) = 0, (27) 

in real space. To obtain this equation from the difference be- 
tween Eq. (IB. 1 21 > and Eq. ( IB. 131 ). we have used rj2q 2 <K b and 
rjq <K b (see below Eq. ([8])), and also \2e q 2 h\ <k \f a - o-\<f>/k\, 
which holds for r » 1 . The last relation can be shown using 
Eq. d24"l) in the long-time limit. Using Eq. ( f27l i. the antisym- 
metric convective mass current j a = j + - j can be expressed 
to first order as 



Ja=Po(v -v ) = -Po^V(r fl - — <py 



(28) 



2(2 2 ^ 3 - l) 1/z zo- These asymptotic behaviors, both for the short- 



vl/2 



this current has a diffusive form. Combining this with the mass 
conservation relations in Eq. (IB. 141) finally yields the diffusion 
equation 

^•-s^-t'H- (29) 

Hence, the local asymmetry in density between the two mono- 
layers finally relaxes by spreading diffusively, with an effective 
diffusion coefficient k/(2b). This slow relaxation is due to in- 
termonolayer friction. Note that we can now interpret the di- 
mensionless number A defined in Eq. (f20b as (half) the ratio 
of the reagent diffusion coefficient D to this effective diffusion 
coefficient. 

The fact that a local density asymmetry in the two monolay- 
ers of a membrane relaxes by spreading diffusively is generic. 
Here, the asymmetry is due to a local chemical modification of 
one monolayer, but it can also be caused, e.g., by a sudden flip 
of some lipids from one monolayer to the other 12 ill . Our de- 
scription of this diffusive behavior generalizes that of Ref. [21], 
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which focused on a perturbation with a spherical cap shape and 
a uniform density asymmetry. Note that this diffusive behavior 
was first mentioned in Ref. B13I1 . within a different theoretical 
framework, and restricting to axially-symmetric deformations. 

3.2. Spontaneous curvature change only 

Let us now focus on the second extreme case, where S] = 0, 
i.e., where only the spontaneous curvature of the upper mono- 
layer is affected by the chemical modification. Fig. [3] shows 
that the deformation converges to a deformed profile, in con- 
trast with the previous case: while a local density asymmetry 
spreads on the whole membrane, a local spontaneous curvature 
modification leads to a locally curved equilibrium shape. The 
deformation is upwards if G > 0, and downwards if G < 0. 
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Figure 3: Dynamics of the membrane deformation in the extreme case Ej = 0, 
where only the spontaneous curvature is changed. The values taken for the other 
dimensionless numbers are those in Eq. d2U . Both the realistic case (i) where 
reagent diffusion is accounted for, and the simpler case (ii) where the chemical 
modification is imposed instantaneously, are considered, (a) Logarithmic plot 
of the height of the membrane deformation H and of the fraction <C of modified 
lipids in front of the reagent source versus dimensionless time t. Both in case 

(i) and in case (ii), H is plotted in units of the extremal value it attains in case 

(ii) . (b) Logarithmic plot of the width W of the membrane deformation and of 
the width of the fraction of modified lipids versus r. Both W and W<j are 
plotted in units of z u . It can be seen on graphs (a) and (b) that cases (i) and (ii) 
yield similar dynamics for r 2> t c ss 0.5. The thin red (gray) lines correspond 
to the analytical laws mentioned in the text. 

As in Sec. 13.11 in the realistic case (i), the membrane dy- 
namics is governed by the reagent diffusion at short times t <k 
t c , while at long times t » t c , it is similar in case (i) and in case 
(ii). The long-time dynamics thus corresponds to the response 
of the membrane to an instantaneously imposed chemical mod- 
ification. 



In the simple case (ii), the approximations on j\ and yo in- 
troduced in Sec. l3.1.2l vield 



h(p,T)K(l-e-"» T/4 )0 s (p), 



(30) 



where we have used the notation p = qzo as above. Calculating 
the inverse Fourier transform of this function for r — yields 



Hoc 



[IT + 4 



(31) 



for all t > 0. This analytical law for the deformation height H 
is plotted in red (gray) lines in Fig. [3} a). 

Let us focus on the long-time limit, which is especially in- 
teresting, since the results found in case (ii) also apply to the 
realistic case (i). Eq. (l30l l shows that in this limit, h(p,r) oc 
4>Ap), so that the long-time profile of the membrane deforma- 
tion is fully determined by that of (p. In particular, for r — > oo, 
W -> 2 V3z and W -> 2 V3z - 

For t — > 0, we find W — > 2(2 2 ^ 3 - l)^ 2 zo again. These 
asymptotic behaviors regarding the deformation width W are 
plotted in red (gray) lines in Fig. [2b), both for the short-time 
limit and for the much more interesting long-time limit. 

As in Sec. 13. 11 in case (ii), p is the only relevant parameter in 
the dynamics of the membrane deformation (see Eq. (l30b ). and 
it determines the transition time between the two asymptotic 
regimes. In the realistic case (i), the value of A is relevant too. 

4. General case 

Generically, a chemical modification will affect both the 
equilibrium density and the spontaneous curvature of the upper 
monolayer. The general solution of Eq. (0 is a linear combi- 
nation of the two solutions obtained for the two extreme cases: 
G — (see Sec. 13. U and 2i = (see Sec. 13.2b . In this Section, 
we will only discuss the realistic case (i) where the reagent dif- 
fusion is taken into account. As in the two extreme cases, the 
effect on the membrane dynamics of the evolution of (p due to 
reagent diffusion is crucial for r <k t c , while it becomes negli- 
gible for r » T c . 

For a lipid membrane in water such that <t > 10~ 8 N/m 
and for zo ^ 5 pm, we have seen that in both extreme cases, 
the dynamics is influenced only by p and A (see Sees. 13. 1|43.21 >. 
Hence, in the general case, the dynamics is influenced by p, A, 
and by the parameter 



KCq 

cr\e 



" L 2 54 



(32) 



which quantifies the relative importance of the spontaneous cur- 
vature change and of the equilibrium density change of the ex- 
ternal monolayer due to the chemical modification. Indeed, 
Eq. (O shows that the ratio of the normal force density aris- 
ing from the spontaneous curvature change to the normal force 
density arising from the equilibrium density change is equal to 
-2a. The effect of varying A was discussed in Sees. 13.1143.21 
and in addition, A cannot vary much for a reagent injected in 
water above a lipid membrane. Hence, we will focus on the 
influence of p and a on the membrane dynamics. 
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The value of the ratio a of the spontaneous curvature change 
to the equilibrium density change is a priori unknown, as <T\ 
and Co, which are involved in Y,\ and G, respectively (see Eq.l20t. 
are unknown. The actual values of cr\ and Co depend on the 
reagent as well as on the membrane itself, as these two param- 
eters describe the linear response of the membrane to a reagent 
(see Eq. 0. Let us assume that a > 0, i.e., that the equilib- 
rium density change and the spontaneous curvature change in- 
duce deformations in the same direction. This is true, e.g., for a 
chemical modification that affects the lipid headgroups in such 
a way that it yields an effective change of the preferred area per 
headgroup [5]. A rough microscopic model, where lipids are 
modeled as cones that favor close-packing, then yields a w 1, 
which corresponds to the case where the two effects yield desta- 
bilizing normal force densities of similar magnitudes [5]. Thus, 
we expect a as 1 . 

Fig-El a ) shows the evolution of H in case (i) for three differ- 
ent values of a. The deformation height features an extremum 
H e , and then a relaxation. This behavior is due to the change of 
the equilibrium density (see Sec. 13. 11 1. Conversely, the nonzero 
asymptotic deformation arises from the change of the spon- 
taneous curvature (see Sec. I3.21 >. The relative importance of H e 
to Hco thus depends on a. This means that studying the dynam- 
ics of the membrane deformation in response to a local chemi- 
cal modification provides information on the ratio of the spon- 
taneous curvature change to the equilibrium density change in- 
duced by this chemical modification. 

In Fig.|5b)-(c), the ratio H e IH^ is plotted versus a for dif- 
ferent values of p. Indeed, as mentioned above, the membrane 
dynamics is determined both by a and by p (at constant A). The 
order of magnitude of H e can be estimated assuming an equi- 
librium density change and/or a spontaneous curvature change 
of a few percent: H e is of order 1 - 10 //m for flaccid mem- 
branes, and smaller than 0.1 /mi if <x > 10 N/m. Hence, we 
choose cto e [10~ 8 , 10" 6 ]N/m, which yields p e [1, 100] for 
Zo = 10//m: such values are taken in Fig.|3b)-(c). 

Our study shows that one can deduce a, and thus the relative 
importance of the spontaneous curvature change to the equilib- 
rium density change due to a chemical modification, from the 
measurement of H e /Ha,. This is very interesting, given that 
such information cannot be deduced from the study of static 
and global membrane modifications rf9l l22tl , Indeed, the equi- 
librium vesicle shapes in the area-difference elasticity model 
are determined by the combined quantity 



Afln 



Aa + ~Cq . 



(33) 



which involves both the equilibrium density and the sponta- 
neous curvature. In this expression, A«o is the dimensionless 
preferred area difference between the two monolayers, which 
is related to the asymmetric equilibrium density change, while 
denotes the dimensionless spontaneous curvature of the bi- 
layer, which is related to the asymmetric spontaneous curvature 
change B22I1 . The parameter £ in Eq. ( 1331 ) is a dimensionless 
number involving the elastic constants of the membrane. Note 
that the vesicle shape variations due to global modifications of 
the vesicle environment are usually interpreted as coming only 
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Figure 4: (a): Logarithmic plot of the height of the membrane deformation H in 
front of the reagent source versus dimensionless time r in the realistic case (i), 
for different values of the parameter a, which quantifies the relative importance 
of the spontaneous curvature change to the equilibrium density change. The 
values taken for the other dimensionless numbers are those in Eq. j2U . In each 
case, H is plotted in units of its asymptotic value in the case a = 1 . The lines 
show the asymptotic values Ha, of H for t — > oo in each case, and H e denotes 
the extremal value of H. (b) and (c): Logarithmic plot of H e /H m versus a in 
case (i), for different values of the dimensionless parameter )i, which affects 
the membrane dynamics as well as a. From down to up: fi = 1; ft = 10 5 = 
3.2; ii — 10; ll — 10 1 5 « 32; u = 10 2 . 



from a change of the spontaneous curvature, under the assump- 
tion that the preferred area per lipid is not modified B9l l23[ I2411 . 

In order to determine the ratio a of the spontaneous cur- 
vature change to the equilibrium density change in a practical 
case, it is necessary to know the value of fi (see Fig.|4fb)-(c)). 
However, as p. does not involve any parameter that depends on 
the reagent (see Eq. (l20li). it is possible to compare the effects 
of different reagents on the same membrane, i.e., their values of 
a, even without knowing the precise value of p.. 



5. Conclusion 

We have analyzed theoretically the spatiotemporal response 
of a membrane submitted to a local heterogeneity, in the real- 
istic case of a local concentration increase of a substance that 
reacts reversibly and instantaneously with the membrane lipid 
headgroups. In general, the dynamics of the membrane defor- 
mation is quite complex, as it involves the evolution of the 
reagent concentration profile due to diffusion in the solution 
above the membrane, simultaneously as the response of the 
membrane. We have shown that, some time after the beginning 
of the reagent concentration increase, the effect of the evolution 
of the reagent concentration becomes negligible. 

Studying this regime enables to extract interesting proper- 
ties of the membrane response. We have shown that a local den- 
sity asymmetry between the two monolayers relaxes by spread- 
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ing diffusively in the whole membrane. Intermonolayer fric- 
tion plays a crucial part in this behavior. In addition, we have 
shown how the relative importance of the spontaneous curva- 
ture change to the equilibrium density change can be extracted 
from the dynamics of the membrane response to the local chem- 
ical modification. This is a significant result since such infor- 
mation cannot be deduced from the study of a static and global 
modification using the area-difference elasticity model. 

Our description provides a theoretical framework for ex- 
periments involving the microinjection of a reagent close to 
biomimetic membranes. In Ref. ||7|], we used the theoretical 
model presented here to analyze experimental results correspond- 
ing to brief microinjections of a basic solution in the regime 
of small deformations, and we obtained good agreement be- 
tween theory and experiment. The results of the present the- 
oretical work show that it would be interesting to conduct ex- 
periments with a continuous injection phase, since the ratio a 
of the spontaneous curvature change to the equilibrium density 
change could then be determined for various membrane com- 
positions and reagents. 

In biomimetic membranes as well as in cells, remarkable 
phenomena occur in the regime of larger deformations: cristae- 
like invaginations fl, tubulation [4], pearling C25ll. budding, 
exo- or endocytosis [ 1], etc. To study such phenomena, it would 
be useful to pursue our study in the nonlinear regime. 

Acknowledgements 

We thank Miglena I. Angelova for motivating us to study 
this subject. We thank Miglena I. Angelova and Nicolas Puff 
for interesting discussions. 



given that the membrane midlayer and the monolayer neutral 
surface are parallel surfaces separated by a distance e. Hence, 
Eq. dA.U can be rewritten as 



This result is independent of the curvature c, contrary to that in 
Eq. ( 1A.1I ). Indeed, by definition, on the neutral surface, curva- 
ture and density are decoupled fT3l . while these two variables 
are coupled on other surfaces. 

Eq. ( IA.3b shows that, due to the chemical modification, the 
scaled equilibrium density on the neutral surface of monolayer 
+ is changed by the amount 

<e q = CqW-<e q (0)=^ (A.4) 

to first order. Besides, Eq. (1A.2I ) indicates that the spontaneous 
curvature of monolayer + is changed by the amount 

<5c e q = C e q(0) - C eq (0) = -C <p (A.5) 

to first order. 

One might wonder why the spontaneous curvature c eq found 
by minimization is not simply -co for <p — (see Eq. ( 1A.2| |). 
This is due to the fact that we work on the membrane midsur- 
face, which is more convenient to study the membrane dynam- 
ics. If the monolayer free energy had been originally written 
using variables and constants defined on the neutral surface, 
then the spontaneous curvature found by minimization would 
correspond exactly to the constant c Q ' which plays the part of co 
when everything is defined on the neutral surface. 



Appendix A. Spontaneous curvature and equilibrium den- 
sity 

The spontaneous curvature and the equilibrium density in 
monolayer + can be obtained by minimizing / + /p + with respect 
to r + and c for a homogeneous monolayer with constant mass. 
First, using the expression of / + in Eq. ©, the minimization 
with respect to r + gives, to first order in e: 



cr (Ti<p 

r„ n — 1 ec . 

eq 2k k 



(A.l) 



Then, the minimization with respect to c yields to first order, 
using Eq. dA. lb : 



Ceq = "CO - C O ~ 



o- e 



(A.2) 



with co = co + 2<T\e/K (as defined in the main text). Note that, 
since we assume that r + = 0(e) and ec - 0(e), we must have 
coe = 0(e) and cro/k = 0(e) for our description to be valid for 
the values of r + and c that minimize / + /p + . This property has 
been used to simplify the results of the minimization. 

The scaled lipid density r+ on the neutral surface of the 
monolayer is related to r + through r* — r + + ec to first order. 
This relation arises from the geometry of parallel surfaces lfl2ll . 



Appendix B. Derivation of the dynamical equations 

Appendix B.l. Hydrodynamics of the fluid above and below 
the membrane 

We wish to determine the velocity field W ± = (w ± , W£) in 
the fluid above (+) and below (-) the membrane. This flow is 
caused by the deformation of the membrane and by the lateral 
flow in the membrane: mathematically, it is determined by the 
boundary conditions corresponding to the continuity of velocity 
at the interface between the fluid and the membrane. 

Given the short lengthscales considered, the dynamics of 
W ± (r, z, t) can be described using Stokes' equation. Adding the 
incompressibility condition, we have: 



j](-q 2 + dfyfy* = (iq + e^P* , 
iq ■ M>* + d z W* = , 



(B.l) 
(B.2) 



where P ± is the excess pressure field in the fluid: the total pres- 
sure is 'P ± = Po + P ± , where Pq is a constant and P ± goes to 
zero far from the membrane. As in the main text, hats indicate 
two-dimensional Fourier transforms. Taking the scalar product 
of q and Eq. ( IB. Il l, and using Eq. ( IB. 2b yields 



P ± = -rjdffi + 4:d 3 W* . 

- q 2 



(B.3) 
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Taking the scalar product of e z and Eq. ( IB. II ). and using Eq. ( IB. 31 ) 
yields 

q A W* - 2 q 2 d 2 z W* + = . (B.4) 

Solving these equations, with the boundary conditions at infin- 
ity: W* — > for z — > ±oo, and at the water-membrane inter- 
face: Wf(q,0,t) = d t h(q,t) and $*(q,0,f) = V\q,t), finally 
yields 

P ± = +2T 1 [q(d l h) + iq-v ± ] e* qz , (B.5) 
= {v* - [iq (d t h) ±qv ± ]z) e Tqz , (B.6) 
Wf = [d,h ± [ 9 (3$) + iq • j>±] z } e* qz . (B.7) 

For our study of the dynamics of the membrane, what is 
needed is the stress exerted by the fluid on the membrane, i.e., 
since we are working at first order in the membrane deforma- 
tion, the stress in z — 0. The viscous stress tensor of the fluid is 
defined by 

T a p = -(J* + Po)S afi + J] (dpW% + d a Wfj , (B.8) 

where a e {x,y,z} and f3 e \x,y,z}- Hence, using Eqs. (IB. 51 
153115771 we obtain 

f +( Z = 0) - f; : (z = 0) = -4 77 q d,h , (B.9) 
f? z (z = 0) = *2t 1 qv ± , (B.10) 

where we have introduced the tangential part Tf z = (T^„ !T*) of 
the stress tensor of the fluid above and below the membrane. 

Appendix B.2. Membrane linear dynamics 

The first dynamical equation we use is a balance of forces 
per unit area acting normally to the membrane. It involves the 
normal elastic force density in the membrane given by Eq. (f5]l 
and the normal viscous stresses exerted by the fluid above and 



below the membrane, which are derived in Appendix B.l (see 
Eq. dB~9l T). It reads: 

-{a-Qq 2 + Tcq^h + keq 2 r a + q 2 (f>-4-7]qd,h = 0, (B.ll) 

where 77 denotes the viscosity of the fluid above and below the 
membrane. 

Besides, as each monolayer is a two-dimensional fluid, we 
write down generalized Stokes equations describing the tangen- 
tial force balance in each monolayer. The first force involved is 
the density of elastic forces given by Eqs. (0 and ©. The sec- 
ond one arises from the viscous stress in the two-dimensional 
flow of lipids. The third one comes from the viscous stress 



exerted by the water, which is derived in Appendix B.l (see 
Eq. ( IB. 10b ). The last force that has to be included is the inter- 
monolayer friction 111311 . We thus obtain: 



! k q (r + - e q 2 h - y- 0j - (772 q 2 + 2rjq)v + 
-b(v + -v") = 0, 
-ikq(r~ + eq 2 hj - (772 q 2 + 2 77 q} v~ 
+ b{v + - v") =0, 



(B.12) 
(B.13) 



where v ± denotes the velocity in monolayer +, while 772 is the 
two-dimensional viscosity of the lipids and b is the intermono- 
layer friction coefficient [13]. 

Finally, we use the conservation of mass in each monolayer 
to first order: 

dtr* + iqv ± = 0. (B.14) 
Combining Eqs. (IB. 1 1MB. 14b , we obtain Eqs. (7]|5). 
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